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Abstract. The simplicial rook graph SR(d, n) is the graph whose vertices are 
the lattice points in the nth dilate of the standard simplex in R d , with two 
vertices adjacent if they differ in exactly two coordinates. We prove that the 
adjacency and Laplacian matrices of SR(3,n) have integral spectrum for ev- 
ery n. The proof proceeds by calculating an explicit eigenbasis. We conjecture 
that SR(d, n) is integral for all d and n, and present evidence in support of this 
conjecture. For n < (j), the evidence indicates that the smallest eigenvalue 
of the adjacency matrix is — n, and that the corresponding eigenspace has di- 
mension given by the Mahonian numbers, which enumerate permutations by 
number of inversions. 



1. Introduction 

Let d and n be nonnegative integers. The simplicial rook graph SR(d, n) is the 
graph with vertices 

V(d, n) := js = (xi, ■ . ■ , x d ) : < Xj < n, ^ Xj = n j 

with two vertices adjacent if they agree in all but two coordinates. This graph has 
N = ( n ^ 1 d ^ L 1 ) vertices and is regular of degree S — (d—l)n. Geometrically, let A d_1 
denote the standard simplex in M. d (i.e., the convex hull of the standard basis vectors 
ei, . . . , e<j) and let nA d_1 denote its n th dilate (i.e., the convex hull of nei, . . . , ne^). 
Then V(d, n) is the set of lattice points in nA d_1 , with two points adjacent if their 
difference is a multiple of — ej for some i,j. Thus the independence number 
of SR(d, n) is the maximum number of nonattacking rooks that can be placed on 
a simplicial chessboard with n + 1 "squares" on each side. Nivasch and Lev [TT] 
and Blackburn, Paterson and Stinson [2] showed independently that for d — 3, this 
independence number is \_(2n + 3)/3j. 

As far as we can tell, the class of simplicial rook graphs has not been studied 
before. For some small values of the parameters, SR(d,n) is a well-known graph: 
SR(2,n) and SR(d,l) are complete of orders n + 1 and d respectively; SR(3, 2) 
is isomorphic to the octahedron; and SR(d, 2) is isomorphic to the Johnson graph 
J(d + 1,2). On the other hand, simplicial rook graphs are not in general strongly 
regular or distance-regular, nor are they line graphs or noncomplete extended p- 
sums (in the sense of [SJ p. 55]). They are also not to be confused with the simplicial 
grid graph, in which two vertices are adjacent only if their difference vector is exactly 
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Figure 1. The graph S7?(3,3). 

— Bj (as opposed to some scalar multiple) nor with the triangular graph T n , which 
is the line graph of K n |3J p.23], [B, §10.1]. 

Let G be a simple graph on vertices [n] = {1, ...,n}. The adjacency matrix 
A = A(G) is the n x n symmetric matrix whose entry is 1 if ij is an edge, 
otherwise. The Laplacian matrix is L — L(G) = D — A, where D is the diagonal 
matrix whose (i, i) entry is the degree of vertex i. The graph G is said to be integral 
(resp. Laplacian integral) if all eigenvalues of A (resp. L) are integers. If G is regular 
of degree 5, then these conditions are equivalent, since every eigenvector of A with 
eigenvalue A is an eigenvector of L with eigenvalue 8 — A. 

We can now state our main theorem. 

Theorem 1.1. For every n>l, the simplicial rook graph SR(3,n) is integral and 
Laplacian integral, with eigenvalues as follows: 

If n — 2m + 1 is odd: 
Eigenvalue of A Eigenvalue of L Multiplicity Eigenvector 



-3 4m + 5 = 2n + 3 ( 2 ™) H a , biC 

-2,-l,...,m-3 3m + 5..., 4m + 4 3 P fc 

m - 1 3m + 3 2 R 

m, . . . , 2m — 1 = n — 2 2m + 3 . . . , 3m + 2 3 

4m + 2 = 2n 1 J 

If n — 2m is even: 

Eigenvalue of A Eigenvalue of L Multiplicity Eigenvector 



-3 4m + 3 = 2n + 3 f™" 1 ) H aAc 

-2,-l,...,m-4 3m + 4, ...,4m + 2 3 P k 

m — 3 3m + 3 2 R 

m- 1, . . . ,2m - 2 = n - 2 2m + 2, . . . , 3m + 1 3 Q fe 

4m = 2n 1 J 



Integrality and Laplacian integrality typically arise from tightly controlled combi- 
natorial structure in special families of graphs, including complete graphs, complete 
bipartite graphs and hypercubes (classical; see, e.g., jT3l §5.6]), Johnson graphs [8], 
Kneser graphs [S] and threshold graphs [TU]. (General references on graph eigen- 
values and related topics include [TJ [31 03 [B].) For simplicial rook graphs, lattice 
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geometry provides this combinatorial structure. To prove Theorem we con- 
struct a basis of r( 2 ) consisting of eigenvectors of A(SR(3,n)), as indicated in 
the tables above. The basis vectors H ai & iC for the largest eigenspace (Prop. 2.6 1 are 
signed characteristic vectors for hexagons centered at lattice points in the interior 
ofnA 3 (see Figure^. The other eigenvectors P fe ,R, Q fc (Props. |2lj| |2~9} |2.10| are 
most easily expressed as certain sums of characteristic vectors of lattice lines. 

Theorem together with Kirchhoff's matrix-tree theorem [5J Lemma 13.2.4] 
implies the following formula for the number of spanning trees of SR(d, n). 

Corollary 1.2. The number of spanning trees of SR(3,n) is 



32(2n 



2)1+2 



3)(v) n' a- 

a— n+2 



3(n + l) 2 (n + 2)(3n + 5) 3 



if n is odd, 



32(2n 



2n+2 



3)( 2 ) n a " 

a— n+2 



if n is even. 

we make the following 



U(n + l)(n + 2) 2 (3n + 4) 3 
Based on experimental evidence gathered using Sage 
conjecture: 

Conjecture 1.3. The graph SR(d,n) is integral for all d and n. 

We discuss the general case in Section [3] The construction of hexagon vectors 
generalizes as follows: for each permutohedron whose vertices are lattice points in 
nA d_1 , its signed characteristic vector is an eigenvector of eigenvalue — ( 2 ) (Propo- 
sition 3.1). This is in fact the smallest eigenvalue of SR(d,n) when n > (f). 



lim 



1. 



Moreover, these eigenvectors are linearly independent and, for fixed d, account for 
"almost all" of the spectrum as n — > oo, in the sense that 

dim (span of permutohedron eigenvectors) 
\V(d,n)\ 

When n < (g), the simplex nA d_1 is too small to contain any lattice permuto- 
hedra. On the other hand, the signed characteristic vectors of partial permutohedra 
(i.e., intersections of lattice permutohedra with SR(d,n)) are eigenvectors with 
eigenvalue —n. Experimental evidence indicates that this is in fact the smallest 
eigenvalue of A(d, n), and that these partial permutohedra form a basis for the 
corresponding eigenspace. Unexpectedly, its dimension appears to be the Maho- 
nian number M(d, n) of permutations in &d with exactly n inversions (sequence 
#A008302 in Sloane [H]). In Section 3.2 we construct a family of eigenvectors by 
placing rooks (ordinary rooks, not simplicial rooks!) on Ferrers boards. 

2. Proof of the Main Theorem 

We begin by reviewing some basic algebraic graph theory; for a general reference, 
see, e.g., [B]. Let G = (V,E) be a simple undirected graph with TV vertices. The 
adjacency matrix A{G) is the N x N matrix whose entry is 1 if vertices i and 
j are adjacent, otherwise. The Laplacian matrix is L{G) — D(G) — A(G), where 
D(G) is the diagonal matrix of vertex degrees. These are both real symmetric 
matrices, so they are diagonalizable, with real eigenvalues, and eigenspaces with 
different eigenvalues are orthogonal |BJ §8.4]. 
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Proposition 2.1. The graph SR(d,n) has ("J^ 1 ) vertices and is regular of de- 
gree [d — l)n. In particular, its adjacency and Laplacian matrices have the same 
eigenvectors. 

Proof. Counting vertices is the classic "stars-and-bars" problem (with n stars and 
d — 1 bars). For each x € V(d, n) and each pair of coordinates i, j, there are X{ + xj 
other vertices that agree with x in all coordinates but i and j. Therefore, the degree 
of x is Ei<i<i<n( x « + Xj) = (d-l) Yh=i Xi = (d- l)n. □ 

The matrices A(d, n) and L(d, n) act on the vector space R N with standard basis 
{eijk : (i,j,k) 6 V(d,n)}. We will sometimes consider the standard basis vectors 
as ordered lexicographically, for the purpose of showing that a collection of vectors 
is linearly independent. 

In the rest of this section, we focus exclusively on the case d = 3, and regard n 
as fixed. We fix N :— ("J 2 )' tne number of vertices of SR(3, n), and abbreviate 
A = A(3,n). 

2.1. Basic linear algebra calculations. Define 



X 

Z fe 



^ E e ^> 

j-\-k—n—i 



B n := {X ! ,Y l ,Z 1 :0< l <n}, 



: {X,.Y..Z,:U- / • n L}. 



e ijki 

i-\-k=n—j 

The vectors Xi,Yj,Zfc are the characteristic vectors of lattice lines in nA 2 ; see 
Figure [2] Note that the symmetric group S3 acts on SR(3,n) (hence on each of 
its eigenspaces) by permuting the coordinates of vertices. 



Lemma 2.2. We have 

n n n 

J = E X > = E Y > = E Z * 



a = E*(x l + Y, 



Proof. The first assertion is immediate. For the second, when we expand the sum 



in terms of the e^^, the coefficient on each e^ is i + j + k = n. 
Proposition 2.3. For every i,j,k, we have 



Ae 



ijk 



Xj + Yj + Zfe — 3e i: ,-fc, 



A3 = 2nJ, 
AX, 



AY,- 



AZi = 



(n — i ■ 
{n — i - 
{n — i - 



2)X, 
2)Y, 
2)Z< 



n — j 

E[x, 

n—i 



z, 



X,+Y, 



□ 

(2.1a) 
(2.1b) 

(2.1c) 
(2.1d) 
(2.1e) 
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Proof. Formula (2.1a) is immediate from the definition of A, and (2.1b) follows 
because SR(3,n) is (2n)-regular. For ( |2.1c[ ), we have 

• IX, Yl Ae 'U-k= Y [Xi+Yj+Z k -3ei dtk ] 

j-\-k—n—i j-\-k—n—i 

= (n-i + l)Xi-3 Y + Y \ Y 3+ Z k] 



j-\-k=n—i 
n—i 



j-\-k=n—i 



(n-i- 2)Xi + Y,i Y j + z j 

3=0 



and (2.1d) and (2.1e) are proved similarly. 



□ 



For future use, we also record (without proof) some elementary summation for- 
mulas. 

Lemma 2.4. The following summations hold: 

n—k—1 n—h—X 

Y [4i-2n]=0, Y [4i-2fc-2-n] = (n-2k- l)(n-2fc-2), 

i=k+l i=k+X 

n-j n-j 

Y [4i-2raj =2(n-j-k)(k-j + l), Y [4i - 2fc — 2 - n] = (n - 2j)(n — k— j). 

i=k+l t=fe+l 

Lemma 2.5. The following summations hold: 

n—k n—k 

Y [4* - 2n] = 0, Y [ 4i - 3n + 2fc - 2] = -(n - 2k + l)(n - 2k + 2), 



i—k 
n-j 



i—k 

n-j 



Y [4« - 2n] = 2(j - fc)(-n + j + k - 1), [4i - 3n + 2fc - 2] = (2j + 2 - 4fc + n)(-n +j + k- 1). 

z— A: i— /c 

Having completed these preliminaries, we now construct the eigenvectors of 
SR(3,n). 

2.2. Hexagon vectors. Let (a,b,c) 6 F(3, n) with a,b,c> 0. The corresponding 
"hexagon vector" is defined as 

H a ,b )C := e a _! ,5, c +l ~ e a,b-l,c+l + e a + l,h-l,c — e a+l,b,c-l + e a,6+l,c-l — e a-l,b+l,c- 

Geometrically, this is the characteristic vector, with alternating signs, of a regular 
lattice hexagon centered at the lattice point (a, b, c) in the interior of nA 2 (see 
Figure [2}. 

Proposition 2.6. XTie vectors {H a ^ c : (a, 6, c) € V(d, n), a,b,c > 0} are linearly 
independent, and each one is an eigenvector of A with eigenvalue —3. 



Proof. The equality AH a i 



-3H a .b jC is straightforward from (2.1a). The lex- 



icographic leading term of H a ,b,c is e a _i.b jC+ i, which is different for each (a, b, c), 
implying linear independence. □ 

Proposition 2.7. Let n > 1 and Ze£ "H„ = {H 0; j, )C : < a, 6, c < n}. Then the 
spaces WhL n and M.B n spanned by T~L n and B n are orthogonal complements in K. . 
In particular, dimK,6„ = (™J 2 ) — (™2 ) = ^n, and the set B' n is a basis for RB n 
(and all linear relations on the Xj, Yj, Z,; are generated by those of Lemma 2.2). 
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Figure 2. (left) The graph SR{3,3). (center) The vector Xi 
and the lattice line it supports, (right) Hi 1 1. 



Proof. The scalar product H tt f,, c • X 4 is clearly zero if the two vectors have disjoint 
supports (i.e., i ^ {a — l,a, a + 1}) and is — 1 + 1 = otherwise (geometrically, 
this corresponds to the statement that any two adjacent vertices in the hexagon 
occur with opposite signs in H aj f,. c ; see Figure |2j. Therefore WH n and RB n are 
orthogonal subspaces of R , and dimMS„ < 3n. For the opposite inequality, we 
induct on n. In the base case n = 1, the vectors Xq, Yj, Zq form a basis of M 3 . For 
larger n, let M n be the matrix with columns X n , Y n , Z n , . . . , X , Yq, Z and rows 
ordered lexicographically, and let M n be M n with the columns reordered as 



Xq^ Y ni Z ni X n . Y n —i, Z n _\^ 



X\, Y , Z . 



For example, 



Af 3 







^3 


z 3 


*3 


Y 2 


z 2 


x 2 


Y 1 


Zi 


X x 


^0 


Zo 


003 


1 





1 























1 





012 


1 














1 





1 














021 


1 











1 











1 











030 


1 


1 





























1 


102 

















1 











1 


1 





111 























1 


1 


1 








120 














1 














1 





1 


201 




















1 





1 





1 





210 




















1 


1 











1 


300 











1 




















1 


1 



If a > 0, then the entries of M n in row (a, b, c) and columns Xi,Y i: Zi equal the 
entries of M„_i in row (a— 1, b, c) and columns Xi_i, Yi, Zi respectively. Hence M n 



has the block form 



u 


* 





M n _i 



where the entries of * are irrelevant and 



1 1 
1 



U 



1 
1 1 



Since rank?/ = 3, it follows by induction that rankM n > rankM n _! + 3 = in. 
Using Lemma [2~2} one can solve for each of X„, Y„, and Z„ as linear combinations 
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of the vectors in B' n . It follows that B' n is a basis, and that the linear relations of 
Lemma 2.2 generate all linear relations on the vectors {Xj, Yj, Z,}. □ 



2.3. Non-Hexagon Eigenvectors. We now determine the other eigenspaces of A. 
The vector J spans an eigenspace of dimension 1; in addition, we will show that 



there is one eigenspace of dimension 2 (Prop. 2.8 1 and two families of eigenspaces of 
dimension 3 (Props. 2.9 and 2.10). Together with the hexagon vectors, these form 
a complete decomposition of R N into eigenspaces of A. Throughout, let a and p 
denote the permutations (1 2 3) and (1 2) (written in cycle notation), respectively, 
so that 

a(X i )=Y i , a(Y j ) = Z j , a(Z k )=X k , p(X l )=Y l , p(Y j ) = X j , p(Z k ) = Z k . 

Proposition 2.8. Let n>l and k = [n/2\. Then 

R := Xfc — Yfe — Xfc + i + Yfc+i 

is a nonzero eigenvector of A with eigenvalue n — k — 3 = (n — 6)/2 if n is even, or 
n — k — 2 = (n — 3)/2 if n is odd. Moreover, the &3-orbit o/R has dimension 2. 



Proof. By §IX<\ . . . pie] , 

AR = (n — fc - 2)(X fc - Y fc ) + ]T [y, 



i—k 



x, 



+ (n-fc-3)(Y fe+1 -X fc+1 ) 



E [ 



X,- 



3=0 3=0 

= {n-k- 2)(X fc - Y fc ) + (Y„_fc - X„_ fc ) + (n - k- 3)(Y fc+1 - X fc+1 ) 



n-k- 2)(Xfc - Y fe ) + (Yfc - X fc ) + (n - k - 3)(Y k+1 - X fc+1 ) if n is even, 

n-k- 2)(Xfc - Yfc) + (Yfc+i - Xfc +1 ) + (n - k - 3)(Y fe+1 - X fc+1 ) if n is odd 

n-k - 3) (X fc - Yfc) + (n - k- 3)(Y k+1 - X fe+1 ) if n is even, 
n- k - 2)(X fe - Yfc) + (n - fc - 2)(Y fe+1 - X fc+1 ) if n is odd 

n — k — 3)R if n is even, 
n — k — 2)R if n is odd 



as desired. The vectors R and c(R) = Yfc — Z k — Y k+ i + Z k+ i are linearly 
independent; on the other hand, p(R) = R and R + c(R) + c 2 (R) = 0, so the 
©3-orbit of R has dimension 2. □ 

Proposition 2.9. For all integers k with < k < L^^J > the vector 

ti— k— 1 

P fc :=-(n-2fc-l)(n-2fc-2)Z n _ fc + ^ [2(t- fc- + (2i-n)(X< + Y< 

i=fe+l 

is a nonzero eigenvector of A with eigenvalue k — 2. Moreover, the & 3 -orbit ofP k 
has dimension 3. 
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Proof. The upper bound on k is equivalent to n — 2k — 2 > 0, so the coefficient of 



Z n _fe in Pfc is nonzero, so P k ^ 0. By (2.1c|. . . (2.1e|, we have 



AP k = -(n - 2k - l)(n - 2k 



2) (^(k-2)Z n _ k + J2 [X,+Y 4 ]^J 



E 

i=fc+l 



•2Z, 



2(/ - /,- 1) I (n-i-2)Zi+^2 [: 

n— i 

+(2i - n) (n - * - 2)(X< + Y<) + £ [x, + Y, 
\ 3=0 

(n -2k- l)(n -2k- 2){k - 2)Z n _ fc - (n - 2k - l)(n - 2k - 2)^ Wi + Y* 

i=0 

n— 1 

+ ^ [(2i-n)(n-i-2)(X i H-Y i )+2(i-fc-l)(ri-i-2)Z j 

i=fc+l 
n— 1 n— i 

+ 51 E [(4i - 2k - 2 - n)(X J + Yj) + (4i - 2n)Z^ 

i=fc+l .7=0 



Interchanging the order of summation in the double sum gives 



A~P k = -(n - 2fc - l)(n — 2fe — 2)(fc — 2)Z n _ 



re 

(n-2k- l)(n - 2fe - 2) J] [x< + Y 4 



n— re— 1 

^ [(2i-n)(n-i-2)(X i + Yi) + 2(i-fc-l)(n-£-2)Z i 
j=fc+i 

fe n-fc-1 

£ [(4i-2fe-2-n)(X i +Y J -) + (4i-2n)Z j 

j=0 i=fc+l 
n— fe— 1 n— j 

Y, E [(4*-2*-2-n)(X i +Y i ) + (4i-2n)Z J -" 

j=fe+l i=fe+l 
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Applying the summation formulas of Lemma |2.4| gives 

k 

AP k = -(n - 2k - l)(n -2k- 2){k - 2)Z n _ fc -(n-2k- l)(n - 2k - 2)^ \Xj + Yj 



n — k—1 

+ [(2i-n)(n-i-2)(X i H-Y i )+2(i-fc-l)(ri-i-2)Z j 

i=k+l 
k 

+ ^ [(n - 2k - l)(n -2k- 2)(X j + Y 3 ) 
j=o 

n — k—1 

+ [( 2 3 - n )( fc + 3 - + Yj) + 2(j - n + k)(j-l- k)Z, 

j=k+i 

(n-2k- l)(n - 2k - 2)(k - 2)Z n _ k 

n — k—1 

+ [(2*-n)(A-2)(X i +Y i )+2(t-fc-l)(fc-2)Z, 



=k+l 



n—k— 1 



= (fc-2) ^-{n-2k-l)(n-2k-2)Z n - k + ^ |(2i - n)(Xj + Y») + 2(i—k— l)Zj 
= (k-2)P k 



i=k+l 



verifying that is an eigenvector, as desired. 

Consider the vectors P k , a(P k ), a 2 (P k ) as elements of the vector space M.B n , 
expanded in terms of the basis B' n (see Prop. 2.7). In these expansions, the basis 
vectors Z„_j., X n _&, Y n _/- occur with nonzero coefficients only in P k , a(P k ), 
cr 2 (Pfe) respectively. This shows that these three vectors are linearly independent. 
On the other hand, p(P k ) = P k , so the 63-orbit of Pfc has dimension 3. □ 



Define 



Q fc := (n-2k+ l)(n - 2k + 2)Z k 



n—k 

J2 [(2j - n)(Xj + Yj) - 2(n - j — k + l)Z 3 



J= k 



(n-2k + \){n - 2k)Z k + (2fc - n)(X fc + Y fe ) 

n— fe 

+ J2 [(2j-n)(X i +Y i )-2(n-j-fc + l)Z 3 



(2.2a) 



(2.2b) 



J= fc+i 

Both of these expressions for Q„ will be useful in what follows. 

Proposition 2.10. For all integers k with < k < [^^^l, the vector Q k is a 
nonzero eigenvector of A with eigenvalue n — k — 2. Moreover, the &3~orbit of Q k 
has dimension 3. 



Proof. The statement is vacuously true if n < 2. By ( |2.2a ), the coefficient of Z 
in Q k is (n — 2k + l)(n — 2k). Provided that k < [ 7JL j^\ , we have n > 2k, so thi 



k 

this 
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coefficient is nonzero, as is the vector Q^. Applying (2.1c). . . (2.1e), we have 



3=0 



AQ k = (n-2k + l)(n-2k+ 2) I (n - k - 2)Z fe + ^ [: 



X,+Y, 



( 2/ - ;/ ) | (n - * - 2)(X< + Y,) + ^ [Xj 



■Y,-+2Z j 



/ n—i 

-2(n — i — k + 1) Un-i- 2)Z* + ^ [x^ + Y 3 



n — fe 



= (n - 2k + l)(n - 2k + 2)(n - fc - 2)Z fe + (n - 2fc + l)(n - 2A; + 2) ^ ... , 



X,-+Y, 



n—k 



[(2i - n)(n - i - 2)(X j; + Y*) - 2(n - i - jfe + l)(n - i - 2)Z 

i—k 

n—k n — i 

Yl, H [( 2 * ~ n )( X J + Y J + 2Z J') - 2(n - i - fc + l)(Xj + Y. 



i=fc j=0 



Interchanging the order of summation in the double sum gives 



n—k 

AQ k = (n -2k+ l)(n - 2k + 2)(n - fe - 2)Z fc + (n - 2fc + l)(n - 2fe + 2) [Xj + Y,- 

3=0 

7i — k 

+ £ [(2j - n)(n - j - 2)(X, + Y,) - 2(n - j - fc + l)(n - i - 2)Z, 

fe — 1 n — k 

+ Y [( 4i ~ 2 ") Z J + ( 4 * - 3n + 2fc - 2)(Xj + Yj) 

n—k n—j 

J — 2 — fe 
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Applying the summation formulas of Lemma |2.5| gives 



i—k 



X,-+Y,- 



AQ k = (n-2k + l){n-2k + 2)(n-k- 2)Z k + (n - 2k + l)(n - 2k + 2) 53 ^ , ... . 

7i— k 

+ [(2j - n)(n - j - 2)(X, + Yj) -2(n- j - k + l)(n- j - 2)Z, 



j = k 
fc-1 



53 [(n - 2fc + l)(n - 2fc + 2)(X j + Yj) 
-ft 



+ 53 [ 2 U - k )(~ n + J+ k - !) z j + (2j + 2 - 4fc + n)(-n + j + k — l)(Xj + Yj 



= (n -2k + l)(n -2k + 2)(n - fc - 2)Z fc 

+ 53 [(» - - 2)(2j - fi)(X,- + Yj) - 2{n -j-k+ l)(n - fc - 2)Zj 

(n - 2fc + l)(n - 2k + 2)Z k + 53 [(2j - n)(Xj + Y 3 -) - 2(n - j - A; + l)Z j 
= (»-fc-2)Q fc 



j=k 



as desired. 

We now show that the 63-orbit of Qj. has dimension 3. Since p(Q&) = Qfc, the 
orbit is spanned by the three vectors Q/c, a(Q k ), cr 2 (Qfc)- We consider two cases: 
fc = and fc > 0. 



First, if k — 0, then the expression (2.2a) for Qo becomes (using Lemma 2.2) 



n 

Qo = (n + l)(n + 2)Z + 53 [( 2 J ~ n )( x : + Y i) ~ 2 ( n ~J + 

3=0 

n n 

= (n + l)(n + 2)Z - 53 [ n ( X J + Y j) + ( 2n + 2 ) Z i + 2 53 + Y J + Z J 

= (n + l)(n + 2)Z - (4n + 2)J + 2nJ = (n 2 + 3n + 2)Z - (2n + 2)3 
= Y (n 2 + n)e tja + 53 {-^n - 2)e l]k . 

i+j—n i,j,k : fc^O 

Accordingly we have 

o-(Qo) = 53 (« 2 + ™) e 0jfc + 53 {-^n - 2)e ijk , 

j-\-k=n i,j->k- i^O 

^ 2 (Qo) = 53 + n ) e iOk + 53 ( _2n ~~ 2 ) e ijk- 

i-J-fc=n iijik: j^O 
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Consider the N x 3 matrix with columns Qo, er(Qo), c 2 (Qo). By the previous 
calculation, the 3x3 minor in rows e„ o ; e 0n0 ; e oon is 



n 
n 



-2n- 2 
n 2 + n 
n + n 



n + n 
-2n -2 

n 2 



II 



II 



= -2(n + lf(n + 2) 2 (n- 1) 



-2n - 2 ; 

which is nonzero (recall that n > 2, otherwise the proposition is vacuously true). 
On the other hand, if < k < [(n— 2)/2j, then (2.2b) expresses Q^, o-(Q/ c ), ct 2 (Qa 



as column vectors in the basis B' n . Let a = 2k — n and b = (n — 2k)(n — 2k 
then the 3x3 minor in rows X^, Zfc is 



i); 



a a b 
aba 
baa 



(2k - nf(n -2k- l)(n - 2/c + 2) 2 
which is nonzero because the assumption k < [(n — 2)/2j implies n > 2k + 2. 



a 



3. SlMPLICIAL ROOK GRAPHS IN ARBITRARY DIMENSION 

We now consider the graph S , i?(d, n) for arbitrary d and n, with adjacency matrix 
A = A(d,n). Recall that SR(d,n) has N :— ( n t_7 1 ) vertices and is regular of 
degree (d — l)n. If two vertices a = (at,..., ad), b = (bt, ■ ■ ■ , bg) € V(d, n) differ 
only in their i th and j th positions (and are therefore adjacent), we write a ~ b. 

Let &d be the symmetric group of order d, and let 21^ C &d be the alternating 
subgroup. Let e be the sign function 



e(a) = 



1 for a G Ql d , 
-I for a <£ 2t d . 



Let Ty- € &d denote the transposition of i and j. Note that &d — 2ld U 2[dT"ij for 
each z, j. 

In analogy to the vectors X, Y, Z used in the d = 3 case, define 

/3: /9~a 



That is, Xq is the characteristic vector of the lattice line through a in direction 
— ej. In part 

indexed by a is 



e; - ej. In particular, if a ~ /?, then X^ j) = xi 4j) . Mor eover, the column of A 



Ae a = -(fje a + Yl ( 3 - 2 ) 

^ ' l<i<j<d 
(i i) 

since e a itself appears in each summand X„ . 

3.1. Permutohedron vectors. We now generalize the construction of hexagon 
vectors to arbitrary dimension. The idea is that for each point p in the interior of 
n /\d-i an j sumc i en tly far away from its boundary, there is a lattice permutohedron 
centered at p, all of whose points are vertices of SR(d,n) (see Figure |3j, and the 
signed characteristic vector of this permutohedron is an eigenvector of A(d,n). 
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Figure 3. A permutohedron vector (n — 6, d = 4). 



Proposition 3.1. Let p,w E R N be vectors such that {p + a(w): a E &d} are 
distinct vertices of SR(d,n). (In particular, the entries of w must all be different.) 
Define 



H„ 



E 



Then H P:t0 is an eigenvector of A with eigenvalue 



Moreover, for a fixed w, 



the collection of all such eigenvectors H p w is linearly independent. 



Proof. By linearity and (3.2), we have 



AH p . w = £ i°) \-\o) e P+°(<<>) + J2 X l:n 

ae&d \ \<i<j<d 

' ffE6i l<i<j<d 

2W + E E £ (-)4S W 



i<i<j<dt?ee d 



H P ,»+ E E [^)xS2 (l0) + -^r;,!X 

l<i<j<d a£$l d 



2 

d ITT 

2 I 



(The summand vanishes because e(cr) = —e(o-Tij) and because changing and Oj 

does not change xi' J '.) For linear independence, it suffices to observe that the 
lexicographic leading term of H p w is e p+ ^, where w denotes the unique increasing 
permutation of w, and that these leading terms are different for different p. □ 
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This result says that we can construct a large eigenspace by fitting many con- 
gruent permutohedra into the dilated simplex. Depending on the parity of d, the 
centers of these permutohedra will be points in 7L d or (Z + \) d ■ 

Let d be a positive integer. The standard offset vector in R d is defined as 

w = w d = ((1 - d)/2, (3 - d)/2, . . . , (d - 3)/2, (d - l)/2) € M d . (3.3) 
Note that w g Z d if d is odd, and w £ (Z + |) d if d is even. 
Proposition 3.2. There are 

2 

d- 1 

distinct vectors p such that H p w is an eigenvector of A(d, n) (and these eigenvectors 



are all linearly independent by Prop. 3.1). 



Proof. First, suppose that d — 2c + 1 is odd. In order to satisfy the conditions of 



Prop. 3.1 it suffices to choose a lattice point p = (ai, . . . , a d ) so that ^2 = fi and 
c < a, < n — c for all z. Subtracting c from each gives a bijection to compositions 
of n — cd with d nonnegative parts and no part greater than n — 2c (that latter 
condition is extraneous for d > 2). The number of these compositions is 

■cd + d-l\ fn- ( rf - 1 ^- 2 ) 



d- 1 J \ d- 1 

Second, suppose that d = 2c is even. Now it suffices to choose a point p = 
(ai + 1/2, . . . } a d + 1/2) G (Z + \) d such that a\ + ■ ■ ■ + a d = n — c and, for 
each i, a* + 1/2 + (1 - d)/2 > and a t + 1/2 + (d - l)/2 < n, that is, i.e., 
c — 1 < a\ < n — c. Subtracting c — 1 from each ai gives a bijection to compositions 
of n — c — d(c — 1) = ri — d(d — l)/2 with d nonnegative parts, none of which can 
be greater than n — d+ 1 (again, the last condition is extraneous). The number of 
these compositions is 

d(d - l)/2 + d - 1\ _ fn - (d ~ 1) 2 (rf ~ 2) ^ 

□ 



d- 1 / V d- 1 



The permutohedron vectors account for "almost all" of the eigenvectors in the 
following sense: if Hd.n C R N be the linear span of the eigenvectors constructed in 
Props. |3.1| and |3.2| then for each fixed d, we have 

(d-l)(d-2) ^ 

dimH ( n ~ ) 

nm TTTTi — rr = um , = 1. (3.4) 

n-S-oo \V(d,n)\ n^oo ("J^ 1 ) 

On the other hand, the combinatorial structure of the remaining eigenvectors is not 
clear. 



The next result is a partial generalization of Proposition 2.7 

(i i ) 

Proposition 3.3. Every H p w is orthogonal to every X„ . 
Proof. By definition we have 



/ a: p+cr(w)~a 
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The index set of this summation admits the fixed-point-free involution (/?, u) <-> 
(t/3,T(t), and e{tu) — — e(cr), so the sum is zero. □ 



Geometrically, Proposition |3.3| says that if a lattice line meets a lattice permu- 
tohedron of the form of Prop. |3.1[ then it does so in exactly two points, whose 
corresponding permutations have opposite signs. 

(i 7*1 

Conjecture 3.4. The vectors span the orthogonal complement ofHd.n- 

This conjecture is equivalent to the statement that every other eigenvector of 
A(d,n) can be written as a linear combination of the X„ . For n < (j), the 

conjecture is that the Xa' 1 '' span all of R^. We have verified this statement com- 
putationally for d — 4 and n < 11, and for d = 5 and n = 7,8,9. Part of the 
difficulty is that it is not clear what subset of the X« ought to form a basis (in 



contrast to the case d — 3, where B' n is a natural choice of basis; see Prop. 2.7) 



3.2. The smallest eigenvalue. For a matrix M with real spectrum, let r(M) de- 
note its smallest eigenvalue, and for a graph H, let t(G) = t(A(G)). The invariant 
t(G) of a graph is important in spectral graph theory; for instance, it is related to 
the independence number [6l Lemma 9.6.2]. 



Proposition 3.5. Suppose that d>l and n > (g). Then r(SR(d,n)) 



Proof. By the construction of Proposition |3.2| there is at least one eigenvector 
with eigenvalue —(2) when n > (if). The following argument that —(2) is m fact 
the smallest eigenvalue was suggested to the authors by Noam Elkies. The edges 
of SR(d, n) in direction (i, j) form a spanning subgraph SR(d, n)i j isomorphic to 

K n+ i + K n + K n -i H + Ki, where + means disjoint union. The eigenvalues of 

K n are n — 1 and — 1, and the spectrum of G + H is the union of the spectra of G 
and H, so r(SR(d,n)ij) = — 1. Since the edge set of SR(d,n) is the disjoint union 
of the edge sets of the SR(d, n)^, we have A(d, n) — J2(i j) A(SR(d 7 n) i: j), and in 
general t(M + N) > r(Af) + t(N), so r(SR{d, n)) > - (f) as desired. □ 

The case n < (2) is more complicated. Experimental evidence indicates that 
the smallest eigenvalue of SR(d, n) is — n, and moreover that the multiplicity of 
this eigenvalue equals the number M(d, n) of permutations in &d with exactly 
n inversions. The numbers M(d, n) are well known in combinatorics as the Ma- 
honian numbers, or as the coefficients of the g-factorial polynomials; see |121 se- 
quence #A008302]. In the rest of this section, we construct M(d,n) linearly inde- 
pendent eigenvectors of eigenvalue — n; however, we do not know how to rule out 
the possibility of additional eigenvectors of equal or smaller eigenvalue 

We review some basics of rook theory; for a general reference, see, e.g., j4j. For a 
sequence of positive integers c = (ci, . . . , q), the skyline board Sky(c) consists of a 
sequence of d columns, with the i th column containing Cj squares. A rook placement 
on Sky(c) consists of a choice of one square in each column. A rook placement is 
proper if all d squares belong to different rows. 

An inversion of a permutation n = (iri, . . . , tt^) G &d is a pair i,j such that 
i < j and 7Ti > ttj. Let &d, n denote the set of permutations of [d] with exactly n 



inversions. 
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Definition 3.6. Let tt G &d,n- The inversion word of tt is a = a(n) = (<2i, . . . , ad), 
where 

a i = #{j G [d] ■ i < j and TTi > ttj}. 
Note that a is a weak composition of n with d parts, hence a vertex of SR{d,n). 
A permutation a € &d,n is 7r- admissible if cr is a proper skyline rook placement on 
Sky(ai + 1, . . . , ad + d); that is, if 

x(a) — a(7r) + w — cr(w) = a(n) + id —a 

is a lattice point in raA d_1 . Note that the coordinates of x(a) sum to n, so admis- 
sibility means that its coordinates are all nonnegative. The set of all 7r-admissible 
permutations is denoted Adm(7r); that is, 

Adm(7r) = {a € &d ■ a-i — <Ji + i > Vi = 1, . . . , d}. 

The corresponding partial permutohedron is 

Parp(7r) = {x{a) : a G Adm(7r)}. 

That is, Parp(7r) is the set of permutations corresponding to lattice points in the 
intersection of nA^" 1 with the standard permutohedron centered at a{ir) + w. The 
partial permutohedron vector is the signed characteristic vector of Parp(7r), that is, 

= ^2 e(cr)e x(o .). 

cr(EParp(7r) 

Example 3.7. Let d = 4 and 7r = 3142 € ©d- Then tt has n = 3 inversions, namely 
12, 14, 34. Its inversion word is accordingly a — (2,0,1,0). The 7r-admissible 
permutations are the proper skyline rook placements on Sky(2+1, 0+2, 1+3, 0+4) = 
Sky(3,2,4,4), namely 1234, 1243, 2134, 2143, 3124, 3142, 3214, 3241 (see Figure|4]). 
The corresponding lattice points x{o~) can be read off from the rook placements by 
counting the number of empty squares above each rook, obtaining respectively 2010, 
2001, 1110, 1101, 0120, 0102, 0030, 0003; these are the neighbors of a in Parp(?r). 
Thus = e 20 io -e 20 oi -emo +enoi -e i2o + eoio2+eoo30-e 0003 ; see Figure[5] 



FIGURE 4. Rook placements on the skyline board Sky(3, 2, 4, 4). 



Theorem 3.8. Let tt G &dn an d A = A(d, n). Then F^ is an eigenvector of A with 
eigenvalue —n. Moreover, for every pair d,n with n < („), the set {F n : n g &d,n} 
is linearly independent. In particular, the dimension of the (—n)-eigenspace of A 
is at least the Mahonian number M(d,n). 

Proof. First, we show that the F T are linearly independent. This follows from the 
observation that the lexicographically leading term of F^ is e a ( ff ), and these terms 
are different for all tt G &d,n- 
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FIGURE 5. The partial permutohedron Parp(3142) in S7?(4, 3). 



Second, let a € Adm(7r). Then the coefficient of e^^) in F,r is e(a) € {1, — 1}. 
We will show that the coefficient of e x ( CT ) in AF^ is —ne(<r), i.e., that 

£(v)^2s{p) = -n, (3.5) 
p 

the sum over all p such that p ~ <r and p € Parp(7r). (Here and subsequently, 
~ denotes adjacency in SR(d,n).) Each such rook placement p is obtained by 
multiplying a by the transposition (i j), that is, by choosing a rook at (i,<7;), 
choosing a second rook at (j,(Tj) with > <7j, and replacing these two rooks with 
rooks in positions (i,<Jj) and (j,<Ji). For each choice of i, there are (a^ + i) — cr, 
possible j's, and X)i(°j + * ™ &i) — n - Moreover, the sign of each such p is opposite 



to that of £7, proving (3.5) 



Third, let y = (j/i, . . . , y d ) € V(d, n) \ Parp(7r). Then the coefficient of e x ^ in 
F,r is 0. We will show that the coefficient of e x / a -i in AF„ is also 0, i.e., that 

where N — {p : x(p) ~ y} (1 Parp(7r). In order to prove this, we will construct a 
sign-reversing involution on N. 

Let a = a(ir) and let b = (b 1: . . . , b d ) = (eti + 1 - y x , a 2 +2 - y 2 , ■ ■ .,a d + d-y d ). 
Note that hi < a% +i for every i; therefore, we can regard b as a rook placement on 
Sky(ai + 1, . . . , a d + d). (It is possible that bi < for one or more i] we will consider 
that case shortly.) To say that y £ F^ is to say that b is not a proper 7r-skyline rook 
placement; on the other hand, we have h = (as would be the case if b were 

proper). Hence the elements of N are the proper 7r-skyline rook skyline placements 
obtained from b by moving one rook up and one other rook down, necessarily by 
the same number of squares. Let b(i t q, j 'I r) denote the rook placement obtained 
by moving the i th rook up to row q and the j rook down to row r. 

We now consider the various possible ways in which b can fail to be proper. 

Case 1: bi < for two or more i. In this case N = 0, because moving only one 
rook up cannot produce a proper 7r-skyline rook placement. 

Case 2: b t < for exactly one i. The other rooks in b cannot all be at different 
heights, because that would imply that h < + (2 H h d) < C^ 1 ) • Therefore, 
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either N = 0, or else bj = bk for some j, k and there are rooks at all heights except 
q and r for some q,r < bj = b^. 

Then b(i j" q, j 4- r) is proper if and only if b(i f q, k J, r) is proper, and likewise 
b(i t r 5 J 4- is proper if and only if b(i f T, k 4- q) is proper. Each of these pairs is 
related by the transposition (j k), so we have the desired sign-reversing involution 
on N. 

Case 3: b{ > 1 for all i. Then the reason that b is not proper must be that 
some row has no rooks and some row has more than one rook. There are several 
subcases: 

Case 3a: For some q ^ r, there are two rooks at height q, no rooks at height 
r, and one rook at every other height. But this is impossible because then ^ 6j = 

Case 3b: There are four or more rooks at height q, or three at height q and two 
or more at height r. In both cases N = 0. 

Case 3c: We have bi = bj = bk] no rooks at heights q or r for some q < r; and 
one rook at every other height. Then 

' b(itr,j±q), b{j^r,iiq), b{k^r,ilq), 
b{itr,klq), b{jtr,klq), f r,j I q). 

For each column of the table above, its two rook placements are related by a 
transposition (e.g., (j k) for the first column) and either both or neither of those 
rook placements are proper (e.g., for the first column, depending on whether or not 
bi < r). Therefore, we have the desired sign-reversing involution on N. 

Case 3d: We have bi — bj = q; b^ — bg = r, and one rook at every other height 
except heights s and t. Now the desired sign-reversing involution on N is toggling 
the rook that gets moved down; for instance, b(J f s, k 4- 1) is proper if and only if 
b(j t s i & 1 1) 1S proper. 



N C 



This completes the proof of (3.6 1, which together with (3.5) completes the proof 



that F,r is an eigenvector of A(d, n) with eigenvalue —n. □ 

Conjecture 3.9. If n < then in fact r(SR(d,n)) = —n, and the dimension of 
the corresponding eigenspace is the Mahonian number M(d,n). 

We have verified this conjecture, using Sage, for all d < 6. It is not clear in 
general how to rule out the possibility of a smaller eigenvalue, or of additional 
(—71) -eigenvectors linearly independent of the F T . 



The proof of Theorem 3.8 implies that every partial permutohedron Parp(7r) 
induces an n-regular subgraph of SR(d,n). Another experimental observation is 
the following: 

Conjecture 3.10. For every n € &d,n, the induced subgraph SR(d, n)|p arp ( 7r ) is 
Laplacian integral. 

We have verified this conjecture, using Sage, for all permutations of length d < 6. 
We do not know what the eigenvalues are, but these graphs are not in general 
strongly regular (as evidenced by the observation that they have more than 3 dis- 
tinct eigenvalues). 

4. Corollaries, alternate methods, and further directions 

4.1. The independence number. The independence number of SR(d,n) can be 
interpreted as the maximum number of nonattacking "rooks" that can be placed on 
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a simplicial chessboard of side length n + 1. By [6l Lemma 9.6.2], the independence 
number a(G) of a ^-regular graph G on N vertices is at most —tN/ (S — t), where 
t is the smallest eigenvalue of A(G). For d = 3 and n > 3, we have r = —3, which 
implies that the independence number a(SR(d,n)) is at most 3(n + 2)(n + l)/(4n + 
6) . This is of course a weaker result (except for a few small values of n) than the 
exact value [_(3n + 3)/2j obtained in [11] and [2]. 

Question 4.1. What is the independence number of SR(d, u)? That is, how many 
nonattacking rooks can be placed on a simplicial chessboard? 

Proposition |3.5| implies the upper bound 

, cp ,. x W d(d+l) fn + d-T 
a(bR(d, n)) < ' 



(2n + d)(d- 1) V d-1 

for n > (2) , but this bound is not sharp (for example, the bound for SR(4, 6) is 
a < 21, but computation indicates that a = 16). 

4.2. Equitable partitions. One approach to determining the spectrum of a graph 
uses the theory of interlacing and equitable partitions [7], chapter 9]. Let X — 
{01, . . . , Ofe} be the set of orbits of vertices of G under the group of automorphisms 
of G. For each two orbits Oi,Oj, define f(i,j) = \N(x) n Oj\ for any x e Oi. 
The choice of x does not matter, so that the function / is well-defined (albeit not 
necessarily symmetric); that is to say, the orbits form an equitable partition of V(G). 
Let P{G) be the k x k square matrix with entries f(i,j)- Then every eigenvalue 
of P is also an eigenvalue of A(G) (H Thm. 9.3.3]. 

When G = SR(n, d), the spectrum of P{G) is typically a proper subset of that 
of A(G). For example, when n = 3 and d — 3, the matrix A(G) has spectrum 
6,1,1,1,0,0,— 2,— 2,— 2,— 3 by Theorem |1.1| but the automorphism group has only 
three orbits, so P(G) is a 3 x 3 matrix and must have a strictly smaller set of 
eigenvalues. In fact its spectrum is 6, 1, —2, which is not a tight interlacing of that 
of A(G) in the sense of Haemers. 

Therefore, these methods may not be sufficient to describe the spectrum of 
SR(n, d) in general. On the other hand, in all cases we have checked computa- 
tionally (d = 4, n < 30; d = 5, n < 25), the matrices P(SR(n,d)) have integral 



spectra, which is consistent with Conjecture 1.3 

Question 4.2. Is SR(d, n) determined up to isomorphism by its spectrum? 
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